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Abstract 

The new model of phonon transmission across the interface between two 
crystals is proposed featured by taking into account the mismatch of crystal 
lattices. It has been found that the mismatch of lattices results in phonon 
scattering at the interface even in the absence of defects. As it has been 
shown, at the normal incidence, longitudinally polarized phonons have much 
larger transmission coefficient than that of transversely polarized phonons, 
excluding the special resonance cases. For the quasi one-dimensional case 
the exact solution has been obtained. 


1 Introduction 

When the heat flows through the boundary between two crystals, temper¬ 
ature at the interface experiences a sharp jump. The proportionality co¬ 
efficient between the heat flux and the temperature jump is known as the 
thermal boundary resistance or the Kapitza resistance. The theory of a 
Kapitza resistance attracts recently serious attention of researchers due to 
its significance for practical use [D El El m. The study of thermal transmis¬ 
sion across the crystal interface, in its own turn, attracts the attention to 
the problem of phonon transmission across the crystal interface since, in the 
majority of cases, the phonons are responsible for the energy transfer across 
the interface. 

Various approximations are used in the theory of thermal boundary resis¬ 
tance to describe the dynamics of the crystal lattice in the boundary region 
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near the interface between two crystals. Two basic approaches are so called 
Acoustic Mismatch Model (AMM) and Diffusive Mismatch Model (DMM). 
In the first approach Eli it is suggested that the dynamics of a lattice and, 
particularly, the interfacial coefficients of the phonon transmission and re¬ 
flection could be calculated using the elasticity theory. However, firstly, this 
approximation is suitable only for a calculation of the Kapitza resistance at 
low temperatures since in this case only low frequency crystal oscillations 
are excited, which are properly well described by the elasticity theory. Sec¬ 
ondly, this approximation doesn’t take into account the phonon scattering at 
the crystal interface. The values of the Kapitza resistance calculated within 
AMM occurred significantly higher than the experimentally determined val¬ 
ues 

Second approach [8] suggests, vice versa, that the interfacial scattering 
is very strong, and phonons incident at the interface “forget” their initial 
direction and uniformly scatter in all directions. The Kapitza resistance 
calculated within this approach is higher than the value calculated within the 
AMM, but still doesn’t give the good agreement with the experimental data 
[7]. Besides, the authors of [7] have studied experimentally the dependence 
of the Kapitza resistance on surface roughness. It turned out that the better 
the surface is polished and, hence, the lower the scattering at it, the higher 
the thermal transport through the surface. But this conclusion contradicts 
the results of the DMM calculations. Thus, the DMM approximation could 
be regarded as refuted. 

Thus, it is necessary to develop more accurate model of the phonon trans¬ 
mission across the interface, which takes into account the atomic structure 
of actual crystals. The simplest model of this kind is the one-dimensional 
chain with the interface. Such a model has been detailed in papers PHD]. 
Though, since an actual crystal is the three-dimensional structure, the one¬ 
dimensional chain permits to elucidate only some qualitative properties of 
the lattice dynamics connected with the atomic structure. The paper m 
describes the three-dimensional model of the lattice dynamics near the in¬ 
terface. However, this model does not take into consideration the lattice 
mismatch of two crystals. 

More complicated models are studied by computer simulation. Thus, in 
paper [12] it was found that both the AMM and DMM do not give the proper 
description of the phonon transmission across the interface. It has been also 
found that the longitudinal phonons have the significantly higher interfacial 
transmittance coefficient than the transverse phonons. The phonon trans- 
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mission across the interface is studied with a Green’s function method in 
[13]. It is also suggested that atoms located at the interface are randomly- 
displaced from the position that they would occupy if they were in volume 
of the crystal. 

The role of anharmonicity of oscillations in the phonon transmission 
across the interface is considered in the paper |l4j . In the most recent study 
[15], dealing with the numerical simulation of the lattice dynamics at the 
crystal interface, it has been found that atoms at the interface oscillate with 
frequency exceeding the maximum possible frequency in a given crystal. 

The presented paper proposes the analytical description of dynamics of 
the crystal lattice near the interface, which takes into account the mutual mis¬ 
match of the lattices. The main assumption is that displacements of atoms 
at the interface are not random, but determined by forces affecting them 
from the side of atoms of the adjacent crystal. Generally the exact solution 
in such a model is not possible. At the same time, the essential qualitative 
effects related to the lattice mismatch can be found. In particular, it is pos¬ 
sible to explain why the transverse phonons have much lower transmission 
coefficients compared to the longitudinal ones. It turns out that the lattice 
mismatch, even in the absence of defects, results in phonon scattering at the 
interface. This scattering is not random, but has a certain structure. 

2 The model 

We consider the interface between two crystals having the structure of a 
simple cubic lattice and contacting by crystal surfaces 

(1, 0, 0) (Fig. 1). The axis x is normal to the interface. In each of the crystals 
we take into account the interaction with the atoms of the first and second 
coordination groups. Such a model of a three-dimensional crystal lattice is 
the most simple and well-studied [16] . 

The lattice constant for the crystal on the left (x is less than zero) and 
on the right (x is larger than zero) is a L ,a R , consequently. The constant of 
the quasi elastic coupling with the atoms of the first and second coordination 
groups for the atoms of the left crystal and for the atoms of the right crystal 
is /3f',/3|', consequently. We neglect the changes of the quasi elastic coupling 
constants and lattice constants near the interface. It is suggested that the 
crystals are infinite, the left crystal occupies a half space x < 0 while the right 
crystal a half space x > 0. It is also assumed that interaction between the 
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Figure 1: Shown are the interface between two crystals (viewed from the 
side) and the atomic bonds considered in the model. The lattice constants 
are denoted as 01 , 02 ; the interfacial atoms are numbered. 

atoms is due to short-range forces, and that with the atoms of the opposite 
crystal interact only the atoms lying at the interface. 

We introduce the following numbering: n x is the numbering of atoms 
along the axis x. For the left crystal the numbering goes from minus infinity 
to zero. Atoms in the interface plane are numbered as zeroes. For the right 
crystal the numbering goes from zero to plus infinity, number zero is for atoms 
in the interface plane. n y , n z is the numbering of atoms in the interface plane, 
it goes from the minus to plus infinity. The bold marking denotes a set of 
indexes n = (n x ,n y ,n z ). 

Atoms located near the interface are influenced by an external poten¬ 
tial produced by atoms from the opposite side of the interface. Due to this 
influence, the crystal is deformed, and the equilibrium position of an atom 
changes relative to its positions in an ideal crystal. Following the theory 
of an ideal crystal lattice without interface, we could expand the potential 
energy in the Tailor series in displacements of atoms from equilibrium po¬ 
sitions and take into account only the minor terms. The Hamiltonian thus 
obtained would not possess the symmetry, and the analytical solution would 
be impossible. Instead, we perform the expansion in displacements of atoms 
from the positions which the nearest-to-interface atoms would occupy if they 
would not interact with atoms of an adjacent crystal. At this approach, the 
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terms of the Tailor series containing the first derivative are not equal to zero 
as it would be the case when the expansion is produced in the neighborhood 
of the true minimum of the potential energy. The derivatives of higher than 
third order, i.e. anharmonic terms, are not taken into consideration. The 
following expressions are valid for atoms on both sides therefore the index 
denoted the side is omitted. 

Thus we obtain: 


U{..,r n ,...) = 


dU 


+ 




d 2 U 


^ dr na n,Q 2^ dr n , a dr i >/3 

n,a ’ n,a: l,p ’ 




(i) 


where r n a - is the displacement of n-th atom along the axis a, a — x,y,z ; 
[/(.., r n ,...) is the potential energy as a function of displacement. 

In the equilibrium position the condition Vn, -£r~ = 0 is fulfilled. We 
substitute the introduced expression for the potential energy into Eq. (1), 
and obtain 


dU 

dr nn 


E 

l,/3 


d 2 U 


dr n ,adr 


ryt = 0. 


( 2 ) 


This system of 3 N equations (N, number of atoms) determines the dis¬ 
placements r[ ] , at zero temperature in the absence of atomic oscillations, or 
O-displacements, that are displacements of atoms from the position about 
which the potential energy is expanded into the Tailor series, to the equilib¬ 
rium position. 

Now, let temperature is not equal to zero and atoms oscillate about the 
equilibrium position. Then the total displacement of an atom consists of 
the O-displacement and its displacement due to thermal oscillations: r n a = 
r n a + u n,a- The expression for the potential energy can be rewritten as 


)+?EE { r n,ad~ u n,a)( r l j p~^ u l,y)- 


n,a! 


^Te drn ’ adri ’P 


According to the Newton second law: 


dU 


dU 


TYlU'Yl.Ct 


-E 


d 2 U 


&V n a j dv n a dv\ ^ 


( r n,a “I" ^n,c 


(3) 

(4) 


By using the definition of the O-displacement (2), we have 

d 2 U 


TflUyt 


E 

i,/J 


dr n , a dri.0 


U r 


(5) 
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Thus, with taking into account only the terms of the Tailor series contain¬ 
ing derivatives of not higher than the second order, the O-displacements are 
completely excluded from the equations of lattice vibrations. This result is 
the generalization of the known property of harmonic oscillator: the constant 
external field does not change its frequency. 

Let us define the interface, at which the average O-displacements are much 
less than the lattice constant, as ideal. Interaction of O-displacements with 
lattice vibrations is revealed only if the terms of Tailor series with third and 
higher derivatives are taken into account. Thus, for an ideal interface the 
interaction of O-displacements with lattice vibrations has the same order of 
magnitude as the interaction of phonons with each other and is revealed with 
taking into account the anharmonicity of vibrations. Hence, we can consider 
this interaction by means of perturbation theory. 

Obviously, if the model of an ideal interface could describe properly the 
real interface, the latter should be sufficiently smooth, and the interaction 
between atoms on opposite sides of the interface essentially weaker than 
interaction between atoms of the same crystal. Atoms on opposite sides of 
the interface should not, in any case, form chemical bonds. 

We continue to consider lattice vibrations in the interfacial region us¬ 
ing the ideal interface approximation and with no regard for anharmonicity. 
Ignoring the terms with derivatives higher than the third, we obtain the 
Hamiltonian function which is invariant with respect to displacements of the 
atoms of the left crystal by an amount equal to the lattice constant a L , and 
also, of the atoms of the right crystal by the amount equal to the lattice 
constant a R , along the axis y or z. So, it is also invariant relative to any 
linear combinations of these displacements na L + ma R of each crystal. The 
symmetry group of this kind determines the specific properties of the lattice 
dynamics at the interface between two crystals. 

3 Equation for the quasi one-dimensional case 

Let’s consider first the more simple case where the phonon falls normally at 
the interface, that is, the wave vector components parallel to the interface 
Qy, Qz = 0. 

We define u R a to be the displacement of the n-th atom lying on the left 
side of the interface along the axis a , and u R , H the displacement of the n'-th 
atom on the right side of the interface along the axis j3. The Newton second 
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law for the 1-th atom of the left crystal lying at the interface is then given by 


mu j iQ = 


_ X ^ d 2 U 

n^l ,/3 


u ii3-^2 


d 2 U 


du l,a du l,/3 


U \,P 


y, D 2 U 
^ du£du2 a 

n',/3 ^ ^ ?P 


U n',/3■ 


( 6 ) 


In the right part of Eq. (6) we separate out the term describing the inter¬ 
action with the atoms of the same crystal as the atom under consideration, 
and the term describing the interaction with the atoms locating on the other 
side of the interface. To do this, we take into account that 


d 2 U 

du i,a du i,p 


d 2 U y^ 

du la du n 4 3 V 


d 2 U 


dutadu*'/ 


(7) 


as follows from the fact that the energy does not change if the atoms of both 
crystals are displaced by an equal distance in the same direction. Substitution 
of Eq. (7) in Eq. (6) gives 


mu Xa = - 


v 92 ‘- 

^ du t d <P 
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n.p 


— U 
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R 
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( 8 ) 


The Erst term in the right part of Eq. (8) describes the interaction of an 
atom with atoms of the same crystal, the second one relates to the interaction 
with atoms lying on the other side of the interface. 

We seek a solution in the form of superposition of the incident, reflected 
and transmitted waves. We do not take into consideration the wave scatter¬ 
ing at the interface so that the wave vector components of the reflected and 
transmitted waves, parallel to the interface, equal zero. In this case, there are 
only three reflected and three transmitted waves with different polarization. 
Let us assume that from the left, at the interface falls the wave of unit am¬ 
plitude and polarization 1. Then, for an atom on the left side (not necessary 
at the interface itself) we obtain 

u n,a = ex P (feat) (exp (- iq[a L n x ) e[ Q + Hi exp ( iq{a L n x ) e[ Q + 

+H 2 exp (iq^a L n x ) e L lct + H 3 exp (iq^a L n x ) e 3a ). (9) 


Here indexes 1, 2, 3 denote polarization, Hi j2i3 amplitudes of the reflected 
waves with different polarization, gi )2)3 are ^-components of the wave vec¬ 
tors. ef 23 stand for the polarization vectors, e a for the component of the 
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polarization vector in direction of the axis a. At the interface n x = 0 , so for 
interfacial atoms each exponent in Eq. (9) becomes unity. 

The wave vector values ^1,2,3 are assumed to be known for the given fre¬ 
quency u, because the disperse relations for the waves in a simple cubic 
lattice are known [T 6 j. Near the interface, it also may be that the values 
3 are imaginary, that is, the oscillation is gradually damping inward the 
crystal. It happens when the frequency of the incident wave is larger than 
the maximum frequency of oscillations in the reflected wave with the cor¬ 
responding polarization [TO]. Further on, we assume that the phonon falls 
at the interface from the left, and so the crystal from the side of which the 
phonon falls, will be called as ’’left” for brevity. 

Similarly, for crystal atoms to the right of the interface we have 

u n', a = ex P { B i ex P (-A/f 1 a R n' x ) ef a + B 2 exp - iq%a R n' x ) e^ Q + 

B 3 exp (- iq*a R n' x ) e $ a ), (10) 

where -81,2,3 are the amplitudes of the transmitted waves. The values gf 2 3 
can be imaginary (see the paper ma and discussion in section 7 ). 

Equations (9) and (10) are the solutions of the Newton equations for 
the atoms lying off the interface, since < 71 , 2,3 and u> satisfy the dispersion 
relations for lattice vibrations, derived without regard for the boundary. The 
problem thus reduces to Ending the Ai, 2 , 3 , -£>1,2,3 which would satisfy the 
Newton equations for the interfacial atoms. For this purpose we substitute 
Eqs.(9) and (10) in Eq.( 8 ) and try to simplify the expression obtained. For 
the term describing the interaction of an atom with atoms of the same crystal, 
we can proceed in a way similar to that in the paper of HU- We divide both 
parts of Eq.( 8 ) by exp (iujt), and then for the first term in the right part we 


have 


d 2 U 


^ .* ()iI ^ - ()iI ^' ' 11 

njtl ,/3 OU l,a OU n,p 


U\ a ) — 


d 2 U 


Y q u l q u l ((1 “ ex P (“*?f ‘ a L n x ))e^ + A,-(l - exp (iqfa L n x ))e ^, 

n^n,/3 n - Q n >/3 j 

d 2 U 


Y *SLl ( x - ex P (-*gfaV))ej9+ 


n^Z,/3 


du l,a du n,f3 

V 4 . v 92 u 

2.^ j 2 Le dul' &U T ' n 

j n/1,/3 *> a n ’^ 


(1 - exp (~iq^a L n x ))e^. 

( 11 ) 


We take into account that the frequency and the polarization vector can be 
expressed in terms of the dynamic matrix [16] 


CJ CjQ — 


1 

m 


y. d 2 U 
du^Mg 


(1 - exp (—iqjd 1 (n x - l x ))e jf s = Y D »pe j p{q j ). 


( 12 ) 

In Eq.(ll) the summation is over the atoms of the same crystal only. To 
the right of the interfacial atom of the left crystal, n x = 0, there are no more 
atoms of that crystal, n x < 0, that can be formally taken into consideration 
in the form of the condition l x < n x . We introduce the notation (see Fig. 2) 


KM 


d 2 U 


- Y 

m 1 7^8 du n,« du tfl 


(1 - exp (—iqfa l)) 


¥n,/3 

l x ^Tlx 


(13) 


which is for the dynamic matrix describing the atom lying at the interface. 
Here we can ignore the fact that n x < 0, since for an interfacial atom it was 
taken into account in the expression l x < n x . So we can operate with Eq.(13) 
formally. 

Denote the difference of the two dynamic matrices (one for an atom in 
the depth of the crystal and the other for an atom at the interface) as 


Datfjisij') D a /3(qj) 


n x 

^a/3 


( ) — — \ &YL 

m 1^9 du n,ad u l,p 
lx 


(1 — exp (- 


-iqfa L l x )). 


(14) 
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Figure 2: Schematically displayed is the interaction between atoms which is 
taken into account in a given dynamic matrix. 


We substitute Eqs.(12), (13) in Eq.(8) and transfer the first term of the 
right part of Eq.(8) to the left. The expression thus obtained can be written 
in new notations in the form 


E Aj E fp + & D af ( qi )e L v = - E 


d 2 U 


>/3 dun, a du*, if) 


K',/3 ~ U n„ 


Now we rearrange the right part of Eq.(15) by using the relation 

d 2 U 


-ft n,a/9 ^ ^ 




(15) 


(16) 


which is the matrix of the interfacial interaction, describing the interaction 
of an atom with atoms lying on the opposite side of the interface, where 
n x = 0 is supposed. Unlike the conventional dynamic matrix, the matrix of 
this kind depends on the number of the interfacial atom, because each atom 
at the interface is differently located relative to atoms on the opposite side 
of the interface. 

Substituting Eqs.(10), (16) into Eq.(15), we come to 


E^E = -E / w4+E<W5-- s >bh 

3 P P 3 

(17) 

Then we sum over all n in Eq.(17) and divide the result by a total number 
of atoms N at the interface on the left. The left part of the expression is 
unchanged for it is independent of the number n of the atom. In the right 
part of the expression we obtain 
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( 18 ) 



n 


This is the averaged matrix describing the interaction of atoms across the 
interface. By performing such an averaging, we do not take into account that 
the atoms nearest to the interface differ in their position relative to atoms on 
the opposite side of the interface, and, hence, interact differently with them. 
Taking account of this difference involves the appearance of scattering, which 
calls for mathematical technique presented below in Sec. 5. In this section 
we proceed ignoring the scattering. 

With the simplifications made above (q VtZ = 0 for all reflected and trans¬ 
mitted waves), all atoms with the same coordinate along the axis x are os¬ 
cillating in phase. In fact, performing such an averaging, we do not consider 
oscillations of an individual atom, but oscillations of the whole of the plane 
near the interface of the crystal. The matrix describes forces acting 
on the crystal plane on the left of the interface from the side of the crystal 
plane lying on the right of the interface. The problem in this case becomes 
quasi-one-dimensional. 

Let the relation 


= 5] ( ®D af (± qj ) + /C e )4 (19) 

P 

be the second dynamical matrix. It expresses the forces acting on atoms 
near the interface in terms of amplitudes of waves differing in polarization 
on the left side of the interface. That to write the analogous expression for 
the forces acting on atoms on the other side of the interface, we introduce 

J aj = K <*P e %- ( 20 ) 

P 

With new notations, Eq.(17) can be rewritten in the form 

Aj + I^B, = M^r. ( 21 ) 

where the summation is supposed to be performed over the repetitive index 
j standing for the wave polarization. 

By writing the Newton second law for atoms on the right of the interface, 
one can rearrange the obtained expression in a similar way. A distinction is 
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in the fact that the averaged matrix of the interfacial interaction is defined 
as 


K r 

1S -/3a 




d 2 U 


N' ^ ^ du R , 


PK ’,p 


= —K. 


N' 


cx.j3'> 


( 22 ) 


where N' is a number of atoms nearest to the interface in the crystal on 
the right. The amounts of atoms near the interface are related as inverse 
squares of the lattice constants. Thus we obtain 


A ^=(S) 2 ^’ ( 23 ) 

which represents the fact that the force acting on the right crystal from 
the side of the left one is equal, in magnitude, to the force acting on the left 
crystal from the side of the right one. Finally, we have 

MitA, + i L ai B, = Mtr 

+ I*A S = 0. (24) 


This is the sought-for system of equations. 


4 Theorem on the interfacial interaction and 
the exact solution in the quasi-one-dimensional 
case 

Let us study the matrix of interfacial interaction, K a ^, using the fact that 
the interaction energy of two crystals is invariant in respect to displacements 
by any linear combination of the lattice vectors. We consider first the general 
case when the lattice parameters are incommensurate, that is, their ratio is 
irrational, a L /a R ^ Q. In this case, linear combinations of the lattice vectors 
and, hence, a set of rearrangements converting the system into itself, is dense 
on the plane. It is reasonable to require that the energy be continuous as 
the function of the relative displacement of the crystals. Since the energy is 
constant on the dense set and, besides, is continuous, it is constant. 

To be specific, we will demonstrate that K y y = 0. We displace the first 
crystal by a small distance Y along the axis y. The energy in a new state will, 
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on the one hand, be equal to the initial one, and on the other, be different 
an amount of the work done by the force of interaction between crystals: 

E = E + f K yy ydy K yy Y 2 /2 — 0 => K y , y — 0. (25) 

Jo 

For the other seven components of the matrix K a y (except for K xx ), the 
proof can be performed in a similar way. Thus, we conclude that, with 
initial assumptions, the only nonzero component of the matrix of interfacial 
interaction is K xx . 

We consider now a case of resonance where a L ja R G Q. It is obvious that 
the above reasoning is inapplicable. Let a L /a R = p/q and p/q < 1. The 
interaction energy of two crystals will be described by a periodical function 
of their relative displacement with a period qa L . We expand it and K a g into 
the Fourier series and substitute the result in Eq.(21). The first, largest term 
of the series turns out to be 


/. - /*’ _ _ TT 

0/31 q 2 (a L ) 21 

We drop the other terms and obtain the following assessment for the case of 
equal potentials of the interatomic interaction: 

K a ,3 ~ 1 /q 2 , a L /a R =p/q 

K af3 = 0, a L /a R £ Q. (27) 

Such a function experiences a discontinuity in all rational points, which is 
impossible from the physical point of view. The paradox can be resolved 
from considerations that this function was obtained ignoring the interaction 
with O-displacements which are large near the resonance. Accounting for 
the anharmonic coupling with O-displacements leads to smearing of spikes in 
rational points. As a result, the actual matrix of the interfacial interaction 
will be given by a smooth envelope of the function (26). 

Hence, we come to the following qualitative result: the atomic oscilla¬ 
tions perpendicular to the plane of the interface interact intensively with the 
atomic oscillations on the opposite side of the interface only in the case when 
a ratio of the lattice constants is close to a rational number with a small 
denominator. Otherwise, this interaction is weak. 

This result can be illustrated by an interesting analogy in the nebular 
mechanics, where the close commensurability of revolution periods of planets, 
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i.e., the proximity of a ratio of periods to a rational number with a small 
denominator, leads to the strong long-period perturbation m- 

Using the predetermined theorem on the matrix of the interfacial inter¬ 
action, one can obtain the exact solution in the quasi-one-dimensional non¬ 
resonance case. We write explicitly the expressions for the elements of the 
matrix ®D a p(qj)\ 



®D,,( qj ) = = ®A«fe) = 0. 


(28) 


The remaining three components are also equal zero since the matrix is sym¬ 
metrical. 

For the case of the normal incidence of the wave at the interface, we have 


f =(i,o, Of 
4 = (0,1, Of 
4 = (0,0, if 


(29) 


Notations: Af, = (3 L . Let a longitudinally polarized wave is incident at the 
interface. Substituting Eqs.(27), (28) in Eq.(18), we get 



(30) 


In addition, one can easily see that all nondiagonal components of the matrix 
M equal zero. Hence, when the incident wave is longitudinally polarized, the 
reflected and transmitted waves have also the longitudinal polarization. The 
system of six linear equations (23) transforms into the simple system of two 
linear equations: 


Mx’i Ai - /3 l B 




( 31 ) 
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Whence it follows 


A p l p r - jgrjW 

1 <i + <i + - f> L f> R 

B p* m r+ -m ri 
1 M*+Mg - 01-0*' 

This solution coincides with the result obtained in the paper [9] for a one¬ 
dimensional chain if one rewrite it in the above notations. In the one¬ 
dimensional case, the 3- by-3 matrices are changed by a single number, and 
f3 L = f3 R . Hence, the results obtained for a one-dimensional chain would be 
applicable in the three-dimensional case. 

Further on we consider the case where the wave with transverse polar¬ 
ization is incident at the interface. For the sake of definiteness, we assume 
that the oscillations occur along the axis y. Since according to the theorem 
on the matrix of interfacial interaction we have K yy = 0, the corresponding 
equations for the transmission coefficient have the form 

My,2 +A 2 = -<f 

M r+ B 2 = 0, (33) 



whence B = 0, \A\ = 1. 

Actually, however, K yy is small but not equals zero identically, because 
the O-displacements, though being small, are nonzero. It was shown in the 
paper [9] that at small frequency of the incident wave, uj —> 0, the trans¬ 
mission coefficient is independent of the interaction force of atoms at the 
interface, but determined by acoustic impedances of media on opposite sides 
of the interface, in other words, the transmission coefficient can be found 
from the elasticity theory. However, when the interaction force of interfacial 
atoms is weak, the transmission coefficient of phonons decreases fast with 
growing frequency. Since the heat transport is carried out by phonons of all 
frequencies, we can neglect the contribution to the energy transport across 
the interface from atomic oscillations parallel to the interface. 
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5 Fourier transform of the matrix of interfa¬ 
cial interaction 

In order to study the phonon transmission across the interface between two 
crystals in more general case when a phonon is incident at arbitrary angle, 
and therewith take into consideration the scattering, the special mathemat¬ 
ical apparatus is required which is presented below. For simplicity we con¬ 
sider a one-dimensional case, becouse the generalization to more dimensions 
is trivial. 

Let we have an infinite one-dimensional chain of atoms spaced a L apart. 
The atoms of the chain are numbered from minus infinity to plus infinity, 
and the axis x is located so that the zero-th atom has the coordinate x = 
0. We place the chain into the external potential $(a;) with a period a R , 
&(x + a R ) = <h(a;). 

Let’s expand the <h(a;) into the Fourier series: 

a R 

$ k =-L- I $(x)e- 27Tikx/aR dx 
0 

+oo 

${x)= $ k e 2nikx/aR 

k=— oo 

Then the potential energy of n-th atom, will equal f(na L ), or 

+oo 

<F„ = e 2 ™ kn ( aL / aR ) 

k=— oo 

One can demonstrate that the inverse rearrangement is also valid: 

1 N 

<L fc = hm — V (36) 

V->oo 2N 

n=—N 



(34) 
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Actually, 


1 


N 


lim —— $ r 

N— 400 2N ‘ ^ 


Z) —2'jrikn(a L /a R ) _ 


n=—N 


= lim —— 

N^-oo 2 N 


N +00 

E(E * 


^2irik'n(a L / a R ) j ^—2 / jrikn(a L / a R ) _ 


n=—N k '=—oo 

+oo 1 N +oo 

V $ fc , lim — V e^V-W-la*) = 

* ^ N—^OO 2N ' ^ 

k '=—oo n=—N k '=—oo 


y ^k'fik'k — (37) 


Thus, we can “forget” the initial function and consider only the 

discrete sets of values <3> n and & k which are expressed in terms of each other. 

For what follows, Eq.(35) can be conveniently rewritten so that it would 
correspond to the expansion in exponents with the wave vectors q k lying in 
the first Brillouin zone, i.e., q k G (—Ti/a L ,n/a L ). To do so, we introduce 


% 


a L \ a R J ’ 


(38) 


where the curly brackets {...} denote the fractional part. Then we have 


$ 


n 


+oo 

y § k e< aLn 

k=—oo 


(39) 


Evidently, this expression is equivalent to Eq.(35), but q' k G (0, 2ir/a L ). That 
to obtain the desired range of values of the wave vector, we have to take 



and, finally, we obtain 


+00 


= y e 

k=— OO 


iqk aL n 


(41) 


Eq.(35) can be applied to the matrix K na p and that to use it for descrip¬ 
tion of the phonon transmission across the interface of two crystal we have 
to rearrange it. Really, 


An,a/3 


d 




dU 


du n,a y du*,f}' 


(42) 
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The function 


dU 


(43) 


4 > = 


n' 


9 K',P '' 


is periodical (the period a R ), since the potential produced by atoms of the 
right crystal for atoms of the left crystal is periodical. Hence, we can intro¬ 
duce 


Ak ’ a/3 N^lo N ^ Kn ’ afj 


0 2ni]<.n.(a L / a R ) 


+oo 


-^-n ,a/3 ^ ^ ^k,a/3 & 


2'Kik.n(a L /a R ) 


k — oo 


(44) 


where k = (0, k y , k z ). 

6 Equation and the theorem on the interfa¬ 
cial interaction in the general case 

Let the wave of the unit amplitude, polarization 1 and wave vector qx falls 
at the interface. Notations: q|| = (0 ,q y ,q z ) are the wave vector components 
parallel to the interface, and ri|| = (0 ,n y ,n z ) is the atom number along the 
axes y, z parallel to the plane of the interface. We seek the solution for the 
left side in the form of superposition of the incident and reflected waves: 


u n,a — ex P fat) (exp ( -iq[ x n x a L + iq[( n||a Z/ )ef Q + 

+ ^ exp (^j,x n x aL + *(q|| + 27rik/a R )n||a L ) e£ ja J , (45) 

J,k 

And for the right side in the form of superposition of the transmitted waves: 

u n',a — ex P fat) ^ Hk'jexp ( iq^ x n' x a R + i(q|| + 27Hk , /a L )n / ||a jR ) e£ ja . 

3, k' 

(46) 
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We substitute Eqs.(45), (46) into Eq.(8). On rearrangement of the left 
part and of the first term in the right part we have 

expf(q||n||a L ) ^exp (2nikn\\a L /a R )A kJ ^ ® D ap{~ q^ x + Q|| + 2irik/a R )e R /3 + 

k,i p 


expi(q||n||a i ) s L) Q/3 (q)ef i/ 3 = 


d 2 U 


(u R , a - Ur 


- du L du R y n ’’p 

t p <JU n,a UU n',0 

(47) 


The Erst term of the right part of Eq.(47) can be written as 

d 2 U 


E 


— dK, a du% 


U n,0 = K n,a0 U n,0, 


(48) 


P 


since u is independent of n' and could be taken out of the sum over n'. 

In the second term of the right part of Eq.(47) we make the following 
rearrangement: 


E 


d 2 U 


n R - 

LL n / a — 


^ du^du^ n 


exp(?;q||n||a L )^ 


d 2 U 


n , 7 p 9 <,a d K ',0 


exp(—iq||n||a i )uf, )/3 . 


(49) 


Let’s divide Eq.(47) by exp i (q||n||a L ) and multiply by exp (27ril/a' R )n||a i ). 
Then we sum over n so that the total number of additives would equal N, 
and divide by N. The left part acquires the form 


'y " -^1-k ,j y ^ ^k,a/3 e k ,j0 E hio y ^ Ko,aP e o t ip- 

kj 0 0 


(50) 


For the second term we introduce 


-^n,a^( c l|t) — y^< 


d 2 U 


n%f} d <« d <,p 


e *qii( n 'ii a ' R ~ n ii ai ) ) 


(51) 




d 2 U 


; - a 9 K a d K> 0 

n'||ni|,p 11 


e *qiiC 


n na —ni 


I a L ) g 2 - 7 rz(k| | ti| | a L /a R + k' 11 n' 11 a H / a L ) 


(52) 
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In these notations the second term of the right part of Eq.(47), after the 
rearrangement, can be written as 

^ -®k ',j -^"lk',a^( < l||) e k , ,j/3- (53) 

3 , k' /? 

Similarly, as it has been made in section 3, we introduce the second 
dynamic matrices: 

the matrix describing the interaction of the waves with the same q| |: 

M k,S = 5Z (+ q|| + 27T?;k/a R ) + K 0 , a p) e£ J/3 , (54) 

P 

the matrix describing the interaction with the oscillations on the opposite 
side of the interface: 


^k,a?( C l|l) — ^ -^kk , ,qg( C l||) e k , ,7'g- (55) 

P 

and the matrix describing the interaction with the oscillations of the same 
crystal, but with another q||: 

Sk ’ai ~ 'y (56) 

P 

In these notations Eq.(47) will be rewritten as 

Mp+Ay, + W + J2 #«(q||)£k'j = <^10 (57) 

k k' 

For the right crystal a similar approach gives: 

+ E + E = 0 (58) 

k' k 

Eqs.(59) and (60) taken together form the complete system of equations 
describing the transmission, reflection and scattering of phonons at the in¬ 
terface between two crystals. 

Thus, it is seen that, on the one hand, even in the model of an ideal 
interface without defects and roughness, there appears the scattering due to 
the mismatch of the crystal lattices. On the other hand, the assumption 
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in DMM that a phonon incident on the interface “forgets” its initial direc¬ 
tion is, generally saying, incorrect. The scattering reveals quite determined 
structure: the wave vector of the transmitted wave, parallel to the interface, 
differs from the wave vector of the incident wave by an integer number of 
the vectors of the reciprocal lattice of the left crystal, while the wave vector, 
parallel to the interface, of the reflected wave differs from the wave vector 
of the incident wave by an integer number of the vectors of the reciprocal 
lattice of the right crystal. That the vectors of the scattered phonons lie in 
the first Brillouin zone it is remained to rearrange them using Eq.(40). 

It is impossible to solve the infinite system of Eqs.(59) and (60) ana¬ 
lytically. However, if the function <3>, in terms of which the matrix K n a p 
is determined in Eq.(42), is rather smooth, the coefficients of its expansion 
into the Fourier series fast decreases and along with them the matrices A'k , a /3 
describing the scattering also do. In this case one can choose the definite 
number of the terms from the system (59, 60) and solve it numerically. In 
the case when the matrices A^ k a/ g are large even at large values of k, one can 
consider the assumption of the DMM as being valid. 

In the case when the scattering is small, one can alternatively use the 
perturbation theory, by taking the solution of the system of six equations 
describing the phonon transmission without scattering as the unperturbed 
one: 


M 0,ajAo,j + -fo,aj(q \\)Bo,j ~ M 0,al 

Mo’ a jB 0,j + IlMoj = 0. (59) 

This system describes the transmission of phonons across the interface with 
no regard for scattering, i.e. the refraction of phonons at the interface of two 
crystals. 

Unfortunately, even this system cannot be solved analytically, since it is 
impossible to inverse the expression (12) and express the wave vector com¬ 
ponent normal to the interface, q x , in terms of frequency and two other 
components (except for the case of high symmetry, q y — q z — 0, consid¬ 
ered in section 4). Such an inversion can be performed for the interface of 
face-centered cubic lattices and was made in the work HU- In the present pa¬ 
per we deal with the model of simple cubic lattices because just this model 
ensured the simplest qualitative study of the impact of the crystal lattice 
mismatch on the transmission of phonons across the interface. 

To facilitate the numerical calculation of the problem, it is reasonable to 
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invoke the theorem on the matrix of the interfacial interaction in the general. 
Let’s expand the function $ into the Fourier series: 

+OO 

$(x)= $ k e 2miix/aR . (60) 

k=—oo 

On substitution into Eq. (42), we find that the matrix K£ a p is expressed in 
terms of the components of <f> as 

^k.a/3 = k a <&k- (61) 

Thus, K(o,k z ), y /3 = 0 and K^ y ,o), z /3 = 0. We use the fact that according to 
Eq.(23), K R a = (a L /a R ) 2 K^p and that for K R a the expression analogous to 
Eq.(42) is valid. We obtain 

Ko^p = f3S ax 5p x . (62) 

Thus, the atomic oscillations parallel to the interface do not contribute to 
the transmission of phonons without scattering at the interface. 

7 Refraction of phonons at the crystal inter¬ 
face 

Eq.(60) admits the analytical solution in the case when a phonon falls at the 
interface at a small angle. Then, by taking q|| as a small parameter and ap¬ 
plying the perturbation theory, we can find the correction to the amplitudes 
of the reflection and transmission, which is proportional to gjj. This solution 
is extraordinary cumbersome and so we omit it here. 

Instead we will give a qualitative description of refraction of phonons at 
the crystal interface, for any angle of incidence. To do so, we display the 
surface of the constant frequency c o in the Brillouin zone (Fig.3) and make a 
projection of this surface onto the plane ( q y ,q z )• If the component q|| of the 
incident phonon is lying in limits of this projection, it is easy to see that for 
the transmitted phonon the wave vector component normal to the interface, 
q x , is equal to the q x of that point in the Brillouin zone which was projected 
into the point qq. 

If the component q|| of the incident phonon lies beyond the projection, 
the transmitted phonon does not appear, because there is no the q x such that 
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Figure 3: A plane of constant frequency in the Brillouin zone. It is shown 
that a phonon with the component of the wave vector parallel to the interface, 
q\\i passes through it. When an incident phonon has the component of the 
wave vector parallel to the interface q\\ 2 , on the opposite side there appears 
an inward-damping oscillation with the damping coefficient k. 

the phonon with the wave vector (q x ,q\\ y ,q\\z) has the frequency u. In this 
case there occurs the oscillation directed into the depth from the interface, 
with the damping coefficient k and wave vector parallel to the interface q||. 
The specific cases of such oscillations are well known Rayleigh waves | [18] . 
Just as there are three branches of oscillations for the given oj, q||, one should 
perform the above procedure for each of them. 

If the frequency of the incident wave exceeds the maximum possible in 
the given branch, there occurs in the crystal the inward-damping oscillation 
at any value of q||. The oscillations of this kind were predicted in the model 
of a one-dimensional chain in paper B33- Probably, it is such oscillations that 
were found in numerical simulations in work [15] . 

The maximum frequency of the transverse oscillations for a simple cubic 
lattice lies inside and not on the boundary of the Brillouin zone p6j. At some 
frequencies, the surface of constant frequency for these oscillation branches 
is closed (Fig. 4). In this case to every q||, lying inside the projection of the 
isofrequency surface onto the plane ( q y ,q z ), correspond two values of q x . To 
these two values of q x correspond different values of the derivative du/dq x , 
and, hence, different directions of propagation of energy. 

The transmitted waves should have the same sign q x as the incident (and 
the reflected the opposite) due to the principle of causality. In addition, 
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Figure 4: It was shown that to each component q||, lying inside the projec¬ 
tion of the isofrequency surface onto the plane (q y , q z ), correspond two values 
of q x . To these two values of q x correspond different values of the derivative 
du/dq x , and, hence, different directions of propagation of energy (shown by 
solid arrows). 

because of the energy conservation law, the energy flux direction is the same 
for the incident and transmitted waves. If we assume that the incident wave 
is one in which the energy propagates towards the interface, then the value 
of q x in such a wave can be both positive and negative. Correspondingly, 
one should take the sign of q x in the transmitted and reflected waves. The 
uncertainty is thus eliminated. 

An interesting effect occurs when the parallel-to-the interface component 
of the wave vector, g||, of the incident phonon is beyond the Brillouin zone 
of the right crystal (Fig.5). In this case the phonon is refracted and under¬ 
goes the Bragg scattering simultaneously. Then it appears that the refracted 
phonons have the opposite direction of propagation, in the plane of the in¬ 
terface, relative to the incident phonons. Thus, the inverse refraction takes 
place. 

It is interesting to note, that the energy flux in the direction perpendic¬ 
ular to the interface should be conserved for each solution of the system of 
Eqs.(60), and for the solutions of the more exact system of Eqs.(58), (59). 
However, for each solution, the energy flux in the plane of the interface can 
be different on the opposite sides from the interface. It means that at the 
interface not only the temperature may undergo a jump but the thermal flux 
parallel to the interface plane as well. 
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Figure 5: Schematically displayed is the surface of the constant frequency in 
the Brillouin zone. It was shown that the incident phonon with the parallel- 
to-interface wave vector component q\ lying outside the Brillouin zone of the 
right crystal, is refracted in the opposite direction. 

8 Conclusion 

In this paper we have studied the model of the interface between two crystals 
taking into consideration the mismatch of crystal lattices. The basic assump¬ 
tion of the model is that displacements of atoms near the interface are not 
random but determined by interaction with atoms of the other crystal. We 
have shown that in the harmonic approximation such the displacements have 
no effect on phonon transmission across the interface. The equation was set 
up, which determines the amplitudes of the transmitted and reflected waves 
of lattice oscillations. The exact solution has been derived for the quasi-one- 
dimensional case. 

It was shown that the mismatch of lattices leads to the scattering of 
phonons at the crystal interface. In other words, the scattering appears even 
in the case of the ideal interface in the absence of defects and roughness. On 
the other hand, such scattering occurs not uniformly in all directions but has 
a certain structure. 

It has been established that the finiteness of the lattice parameter results, 
at certain angles of incidence, in the backward phonon refraction at the 
interface. A new family of lattice oscillations has been described, which 
includes only the atoms lying near the crystal interface. It was predicted that 
at the plane parallel to the interface, the heat flows suffer a discontinuity at 
the interface. 

The main result of our study is that we have shown that the oscillations 
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of atoms in the plane of the interface interact weakly with the oscillations of 
atoms on the opposite side of the interface, except of specific resonance cases. 
In the case of normal incidence of phonons at the interface, the transmission 
coefficient of the transversely polarized phonons is much less as compared 
to that of the longitudinally polarized phonons. For an arbitrary angle of 
incidence, the transmission coefficients of phonons of any polarization are 
less than that calculated with the elasticity theory, even in the case of low 
frequencies. Allowance for this factor leads to that the calculated value of 
Kapitza resistance is approximately three times greater. Calculation of the 
interface thermal resistance, performed by means of the method proposed in 
the work [19], gives significantly lower values as compared to the experiment. 
Accounting for the smallness of the transmission coefficient of transversely 
polarized phonons explains this discrepancy. 

The author thanks E. D. Eidelman and G. V. Budkin for helpful criti¬ 
cal comments, and A. Ya. Vul’ for attention to this study. The work was 
supported by Dynasty Foundation. 
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